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Abstract 

An improved analysis of the b — > s + 7 decay in the minimal flavor violating case is 
given taking into account additional contributions in the supersymmetric sector which 
enter in the next-to-leading-order (NLO) and are enhanced by tan f3 factors. Specifically, 
we compute a set of twenty one-loop diagrams to give the most complete analysis to 
date of the NLO supersymmetric corrections. These modifications are computed from the 
effective charged Higgs and neutral Higgs couplings involving twelve loop diagrams for 
the charged Higgs sector and eight loop diagrams for the neutral Higgs sector. While the 
computations of these corrections are available in the literature, their full forms including 
the complex phase dependence has not be considered. Our analysis takes account of the 
full allowed set of twenty one-loop diagrams and is more general since it also includes the 
full dependence on CP phases in non universal sugra and MSSM models. A numerical 
analysis is carried out to estimate the size of the corrections to b — > s + 7. We also briefly 
discuss the implications of these results for the search for supersymmetry. 



* Current address of T.I. 



1 Introduction 



One of the most severe phenomenological constraints on supersymmetric (SUSY) models 
arises from the measurement of the inclusive rare decay B —> X s ^. This decay only occurs 
at the one-loop level in the standard model(SM)pQ, and therefore the supersymmetric 
radiative corrections are important and might even be of the same order of magnitude 
as the SM contribution (For early work on supersymmetric contributions to b — > 37 
and implications see Refs.j2lE])- in this paper we carry out an improved analysis of the 
branching ratio BR(6 — > 57) assuming the extended minimal-flavor- violation (EMFV). By 
EMFV we mean that the squark and quark mass matrices are diagonalized with the same 
unitary transformation, in which case the only source of flavor violation is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix but CP violation can arise in our model from both the 
CKM matrix and also from the soft susy breaking parameters. The strong constraints on 
flavor changing neutral current, indeed suggest a kind of organizing principle like EMFV 
for the case of softly broken supersymmetry. 

The new results presented in this paper consist of the complete calculation of the 
supersymmetric one-loop corrections to the Higgs sector couplings that enter into the 
calculation of the next-to-leading-order contributions to BR(5 — * 57) through corrections 
to vertex factors. These beyond-leading-order SUSY corrections are parameterized by 
three e's; e&(i), e*(s) and e&& and can have large effects due to contributions that are 
enhanced by factors of tan/3. In this paper we derive the tan/3 enhanced as well as the 
tan /3 non-enhanced contributions. Of course there exist two-loop (NLO) supersymmetric 
corrections beyond the ones parametrized by the e's. However, such NLO corrections are 
generally small or can be absorbed in a redefinition of the SUSY parameters |U E] • As is 
well known the precision theoretical analyses of sparticle masses and couplings are strongly 
affected by the b — > 57 constraint and such predictions would be tested at colliders in the 
future. The above provides the motivation for an improved b — > 57 analysis which is the 
purpose of this analysis. 

The current average value for the BR(6 — > 57) from the experimental data jH] is, 

BR(b — > 87) = (355 ± 24+i ± 3) x 10~ 6 , (1) 

by the Heavy Flavor Averaging Group [Zj. 

The standard model result depends sensitively on the QCD corrections jH] and we will 
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use the value 

BR(B -> X s7 ) = (3.73 ± .30) x KT 4 , (2) 

which takes into account NLO QCD corrections. In this analysis we largely follow 
the analysis of the micrOMEGAs group [TO], i n the computation of the BR(6 — > 57), with 
exception of the calculation of the beyond-leading order SUSY corrections. Further, we 
extend to the case of non-zero CP-phases. In the following we give the essential basics 
of the analysis and refer the reader to the previous literature for more details (see, e.g, 
Ref. PU] and references therein.). The theoretical analysis of b — > s 7 decay is based on 
the following effective Hamiltonian 

AC 8 

H eff = --£KVa y E c i(Q)°i(Q) ( 3 ) 

where Vtb and Vu are elements of the CKM matrix, Oi(Q) are the operators defined 
below and Ci(Q) are the Wilson coefficients evaluated at the scale Q. The only Wilson 
coefficients that contribute are C 2 , C 7 and C 8 and the corresponding operators are defined 
as follows (see e.g., Ref.jH]) 

2 = (c L ^b L )(s L ^c L ) 

7 = -^— 2 m b {s L a^b R )F^ (4) 



Os = ^m b (s L a^T a b R )G% 



Here e is the magnitude of the electronic charge, g s is the strong coupling constant, T a 
(a=l,..,8) are the generators of SU{2>)c and G a ^ v are the gluonic field strengths. As is 
well known the decay width T(B Xg'y) has an ml dependence and thus subject to 
significant uncertainty arising from the uncertainty in the b quark mass measurement. 
However, the semileptonic decay width T(B — > X e eu) also has the same ml dependence 
but is experimentally well determined. For this reason one considers the ratio of the two 
decay widths where the strong dependence cancels out. The ratio of interest including 
the photon detection threshold is defined by jHl ITT] 



F(B — > ^s7)l_E; 7 >(i_5)£;max Q a 



v t * s v tb 



V, 



do 
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K NLO (S) , (5) 



Y{B^X c eV) nf(z) 

where f(z) = 1 — 8z 2 + 8z 6 — z 8 — 24z 4 \nz is a phase space factor and z = {m c /m b ) is given 
in terms of pole masses. We take 5, which is related to the photon detection threshold, 
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to be 0.9 and T(B — > X c eu) to be 0.1045. K^lo depend on the Wilson coefficients and 
is given in the form ^TJ ^] 



K NLO (S) = £ k^Q b )ReC^(Q b )C 0) *(Q b ) + S(5)^^Re[Cf (Q 6 )<^ UJ *(Q 6 )] 

,,,=2,7,8 L 27r 

+ ^)^7T (2Re[C} cm) (g fe )Cf } *(Q 6 )] - A; (cm to)|Cf (Q fe )| 2 ) (6) 
a s {Q b ) v 7 

where kij,S(5) are as defined in Ref.^2]; and we use the running charm mass m c (m b ) 
as suggested in Ref.p. We take the renormalization scale, Q b , to be the b-quark mass. 
Above the Wilson coefficients have been expanded in terms of leading-order and next-to- 
leading order as follows [TT] 

Q(Q b ) = cf\Q b ) + ^f^c\ l \Q b ) + -^ct m \Q b ) . (7) 

The coefficients to leading order at the scale of the b-quark mass can be obtained from 
the Wilson coefficients at the electroweak scale Qw by renormalization group evolution 
such that 



ci°\Q b ) = Uv-l 



_6_ 

4-7^23 



Cf(Q 6 ) = V ^C^(M w ) + ^{ v -s- v -^C^(M w )+J2h iV ^ 



cf(Q 6 ) = ^(C^M^ + ^J+^h^, (8) 

where i] = a s (Mw) / ot s (Q b ) and hi, hi, a, and hi are numerical coefficients and are listed 
in Appendix AjH]. The next-to-leading order contributions and k cm are defined as in 
Ref. [11110] • 

The main focus of this paper is the next-to-leading-order supersymmetric contributions 
to the Wilson coefficients C^s at the electroweak scale. Here Cj$ are sums of the Standard 
Model contribution arising from the exchange of the W and from the exchange of the 
charged Higgs and the charginos, so that 

C7${Qw) = C^ 8 (Qw) + C?s (Qw) + C^8 (Qw) ■ (9) 



Additonally the gluino exchange contribution has been computed in Ref. [T3j. However, 
contributions to the Wilson coefficients arising from gluino and neutralino exchange are 
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negligible in the MFV scenario. Studies of BR(6 — * sj) beyond the MFV scenario, by 
looking at the effects from generational squark mixing, has recently been performed in 
Ref. ^3]. In the analysis of the supersymmetric contributions to the next-to-leading- 
order we will take into account the CP phase dependence. It is now well known that 
large CP phases can appear in SUSY, string and brane models while still allowing for the 
possibility of electric dipole moments of the electron, of the neutron and of the 199 Hg atom 
consistent with experiment |151 IT^l IT71 ITH] . (For the current experiment on the EDMs see 
Refs. |19( I2U1 I2T].) If phases are large they will have important effects on a number of 
phenomena 123 12H 123 123 123 123 123 EOJ - 

The outline of the rest of the paper is as follows: In Sec. 2 we give the effective La- 
grangian for the charged Higgs and exhibit how the corrections e' b (t), e' t (s) and e bb , which 
bring in tan f3 factors, enter in the charged Higgs Yukawa couplings. In Sec. 3 we exhibit 
the dependence on e' b (t), e' t (s) and of the Wilson coefficients Cj$. In Sec. 4 we give a 
comparison of our work with previous ones. A numerical analysis is given in Sec. 5 and 
we determine regions of the parameter space where sizeable differences occur using the 
full formulae derived in this paper relative to the partial results of some of the previous 
works. Conclusions are given in Sec. 6. In appendix A, the parameters hi, hi, and 
hi that appear in Eq. (JBJ) before are listed. In appendix B we give an analysis of e' b (t) 
by computing the six diagrams in Fig.jTJ. In appendix C we give an analysis of e£(s) 
by computing the six diagrams in Fig.(j2J). In appendix D we give an analysis of e&& by 
computing the six diagrams in Fig.Q and the two diagrams of Fig.fjlJ). 

2 Effective Lagrangian 

To discuss the beyond-leading-order supersymmetric contribution it is convenient to look 
at the effective Lagrangian describing the interactions of quarks with the charged Higgs 
fields H ± and the charged Goldstones G ± . We use the framework of the minimal su- 
persymmetric standard model (MSSM) which contains two isodoublets of Higgs bosons. 
Thus for the Higgs sector we have 



The components of Hi and H 2 interact with the quarks at the tree level through 




(10) 



- £ = eijhtbjiHlQ^ - 



eijhttnHiQi + H.c 



(11) 
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The SUSY QCD and the SUSY Electroweak loop corrections produce shifts in these 
couplings and generate new ones as follows 



- C cS = €ij(h b + 5hl)b R HiQl + Ahib R Hi*Qi 

+ 5hi)i R HlQ j L + Ahit R Hi*Qi + H.c (12) 

where the complex conjugate is needed to get a gauge invariant £ e ff- We note that in 
the approximation 

8h) = 5h 2 f , 

Ah) = Ah} (13) 

one finds that Eq. ()12j) preserves weak isospin. This is the approximation that is often used 
in the literature. However, in general the equalities of Eq. (|13|) will not hold and there will 
be violations of weak isospin. It has been demonstrated that the weak isospin violation 
can be quite significant, i.e, as mush as 40 — 50% or more of the total loop correction to 
the Yukawa coupling [29] . 

The typical supersymmetric loop that contributes to the shifts in the couplings Ah l j- 
and 5h l j contains one heavy fermion / and two heavy scalars Si and S2 or one heavy 
scalar S and two heavy fermions fi and ^2- The basic integral that enters in the first case 
is 

j = f dk m f+ # (u) 

1 J (2vr) 4 (k 2 - m)){k 2 - m 2 Si ){k 2 - m%) 1 ' 

The basic integral that enters in the second case is 

d A k {m fl + ^)(m /2 + #) 



(15) 



(2tt) 4 (k 2 - m 2 h )(k 2 - m 2 2 )(k 2 - m 2 ) 

The largest finite parts of these integrals that contribute to the vertex corrections, in the 
zero external momentum analysis would read 

1 22 

11 = T^i — H \ — 2~> — r) 



1 n 9 2 V 2 ' 

\bTi z m s m s 



mt 



where the function H is given by 



(l-x)(x-y) (l-vKv-x) 
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in case x 7^ y and 

H(x,y)=H(x)= - [l-s + lns] (18) 
[x I) 

for the case x = y. 

Electroweak symmetry is broken spontaneously by giving vacuum expectation value 
v x l\pl to H\ and v 2 /V2 to Iff. Then the mass terms for the quarks arising from Eq.([12|) 
would be 

- C m = m b b R b L + m t t R t L + H.c (19) 
with m;, and m t related to h b and h t as follows 



h 



\[2m b 



b 



vi(l + e bb tan(3) 



= ^ (20) 

w 2 (l + e tt cot/?) v ' 



where 



Ahl Jh\ 
e bb = -r 1 + cot/5— ± 

e« = — ^+tan/3— ^ (21) 
h t ht 

The electroweak eigenstates of charged Higgs interaction with quarks in Eq. (jl2|) is 
C eS = {h* b + 8hf)i L b R Hl* - Ahl*t L b R Hl 

+(ht + 5h])t R b L Hl - Ah 2 t i R b L H^ + H.c (22) 
By going from the electroweak eigenstates basis to the mass eigenstate H + and G + basis 

Hi* = sin(3H + - cos (3G + 

H\ = cos f3H + + sin (3G + (23) 

and by using Eq. (J2U|) one gets 
r 5 n+5 l + et(fo)cot/3 l + 4(t)tan/3 

^ff = 7=— G + {m f — Tir t d } L - m b — ^-t L b R \ (24) 

V2M W l + e^cot/3 l + e£ b tan/3 

rr , r l + 4(&) tan 5 _- , l + e 6 (t)cot/3 „- , , 

+ — JL- ff+{ mf ; -f cot (3t R bi, + m b " ; ^ tan/ft^} + if.c. 

v 2Mw 1 + e« cot /3 1 + e£ b tan /? 
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with 

e t (b) = — ^+tan/3— ^ 
rk lit 

Ah 1 * fih 2 * 
n b n b 

Ah 1 * fih 2 * 

e b (t) = -^ + tan/3^- (25) 

For flavor mixing to be considered we should have worked out the analysis for three 
generations of quarks from the beginning. Thus the general effective largrangian would 
read 

r 9 n +f^ T/ 1 + e t (d) cot - ^ 1 + e' b (u) tan (3 _ 

V2M W ^ l + e it cot/5 V l + e£ 6 tan/3 

Q rr j. r ^ T . 1 + ei(<i) tan/3 _- , 
+ —J—H + {£m t V td - , tV ; -f cot pt R d L 26 
v^Mvi/ Y l + e tt cot/3 

^ T ^ l + eb(w)cot/3 , TT 

+ 2^ -— tt tan fJu L b R } + if.c. 

where V qq i here are the radiatively corrected CKM matrix elements. 2 

The terms with e u can be ignored since the radiative corrections for the top quark 
mass are typically less than 1% [21]. As in (3JEI1III2] we will ignore the terms with e t (d) 
and €b{u). This is a good approximation in the large tan/3 region, and for small values of 
tan /3 these e's have little or no influence on the rate of b — > sj. 



3 Wilson Coefficients 

Using the above Lagrangian for the interactions of quarks with the charged Higgs and 
the charged Goldstone bosons along with the Lagrangian that describes the interaction 
of quarks with W ± bosons: 

C = g V td t L 'fd L W+ + H.c, (27) 

d 

2 For a precise analysis of introducing CKM matrix elements into the Lagrangian and their radiative 
corrections, see |5T] , 
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the contributions to C7 g from the iy-boson and from the charged Higgs are given by: 



(%{Qw) = FgM + ^^g^ fffe.) (28) 



where and ?/t are defined by 



m 2 t (Q w ) rn 2 t {Q 



^ = '^ L1 , m='^P- (so) 



and F$ and F$ are given by 

. x(7-5x-8x 2 ) x 2 (3x-2) n 
= 24(,-l)3 

(2) = ^(3 - 5x) x(3x - 2) 

7 v ; 12(x- l) 3 6(z- l) 3 

, x(2 + 5x — x 2 ) 3x 2 

8 V ; 8(x-l) 3 4(x-l) 4 
(2) _ x(3-x) x 

^ 8 ^ " 4(x-l) 3 2(x-l) 3lnX [6i) 

In the limit where all the supersymmetric particles becomes heavy, the SUSY correction 
to the W contribution vanishes. Thus, in this decoupling limit one finds e* hh = e' b (t). 

The chargino exchange contribution to with the beyond-leading-order SUSY cor- 
rections, has been derived in Ref.|3] and extended to the case of non-zero CP-phases in 
Ref.|32|- We have 

k=lt=l t k x f Xi xf 

+ tlhu^F^) + rlf^F^)) (32) 

where Q s is the soft SUSY scale and m\i is the mass of the first and second generation 
up-type squarks, which we take to be identical. Further, 



r . - D* \/* - m *(^ s ) n* V* r' - D tli U j2 

v2Miysmp V2cosp(l + e^tanp) 
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and where and r\- are obtained from and r\- by setting the matrix D t to unity. 
Finally the loop functions F 7 8 (x) appearing in Eq. (j32|) are given by 

The value of the chargino contribution at the scale Qiy is computed as in Ref.jl], where 
we use (3q = —7 corresponding to six flavors. Only the chargino contribution may give a 
CP-violating contribution at the leading order. However, as the e's may be complex all 
three contributions; the W, the charged Higgs as well as the chargino, may be complex at 
NLO order. We note that all the NLO SUSY corrections scales with 1/(1 + e* hh tan f3). To 
complete the analysis what remains to be done is the computation of e' b (t), e' t (s) and e^ and 
as mentioned above we will compute these in the zero external momentum analysis [35] . 
However, we will calculate all one-loop SUSY QCD and SUSY electroweak corrections to 
these for any tan/5. We collect the expressions for these corrections in Appendices B, C 
and D. While analyses for these exist in the literature they are not fully general when CP 
phases are present and the soft parameters are in general complex. 



4 Comparison with previous works 

In this section we compare our results with previous works 3 . We start by comparing 
our results with the work of [23] (BCRS), as this analysis is the most complete of the 
previous works. It includes the exact one-loop results for the e's, but only in the limit 
of CP conservation. Thus BCRS considered the one-loop corrections to the vertex of the 
charged Higgs and charged Goldstones with quarks in their Figure 4. The (AF£) JI and 
(AF*) JI in Eqs. (A.8) and (A.9) of Ref.[M|, where J = 1,2,3 for u,c,t and I = 1,2,3 
for d, s, b are related to our e's as follows: 

(AFl) JI = Vjjh j sin (3e'j(I) 
(AF 2 L ) JI = V JI h J cos (3ej(I) 
(AF^) JI = V JI h* I cos /3ej(J) 

(AF 2 R ) JI = -Vjjh*j sin Pe'j(J) (35) 

3 A brief comparison of partial analysis of the e's with previously works was given in Ref. |28| 
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The first and third lines of the above set are for the charged Higgs H + and the second and 
fourth are for the charged Goldstone G + . The relations relevant for the current analysis 
are the first and the fourth ones and thus we will explicitely check the validity of these. 

In order to compare the e's of our appendices B and C and vertex corrections presented 
in the appendix A. 3 in Ref . |34| . we first establish a dictionary connecting the notation 
in the two works. Thus the form factors in the two works are related to each other by 
xCq(x, y, z) = —H(-, -), where H(^, -) is the form factor used in our work and Cq(x, y, z) 
is the form factor used by BCRS. The diagonalizing matrices Z % 1{Z % +) of BCRS correspond 
to our UjiiVj*) and Z 1 ^ of BCRS corresponds to our Xy. In our analysis we did not consider 
flavor mixing in the squark sector, so the squark mass-squared matrices are 2x2 and 
not 6x6 ones. Thus in the case of e' b (t), the sum in the expressions for AF is over the 
squark mass eigenstates of the third generation; In another words we sum over the third 
and sixth entries in their matrices. So the Zjf* of BCRS corresponds to our D^j, Z[/ 
of BCRS corresponds to our D ti j with I, J — 3i,3j. In the case of e' t (s) we sum, for 
the s squark over the second and fifth entries in their matrices and the Zj/* of BCRS 
corresponds to our D S ij with I, J = 3i — 1, 3j — 1. As BCRS follows the conventions of 
Ref. jni], the Higgs coupling hb in their Lagrangian is our —tit, and their h t is equal to 
ours. The trilinear coupling A t in their superpotential is our —h t A t and their A\, is our 
hbAb as can be seen by comparing the superpotential in Sec. 3 of their reference [30] and 
the superpotential we are using which is the same as in Eq. (4.15) of Gunion and Haber 
|36j . Also the elements Z 1 ^ are defined in section 4 of Ref. |3*o] . 

We now give details of the comparison. The first term in (AFl) JI of Eq.(A.8) of 
Ref . |34j has the following correspondence in our notation 

(AFi)S term - Vjjhj sin (3(e'i l) (I) +ef ] (I)) (36) 

In comparing our results with theirs one finds that we have an explicit gluino phase 
dependence £3 in our analysis. This gives us the maximum freedom in the choice of the 
independent set of phases to carry out the analysis in. 

The second term in (AFl) JI of Eq. (A. 8) of Ref. [3*4] has the following correspondence in 
our notation 

(A^) 2 J nd term - VjshjsmPie'fil) + ef\l)) (37) 
Using e = gsm9 w and g\ = gta.n9w, one can prove that 

(ajkDjtj - -f Jk D J2j ) = — ^V^/tf* 
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{(3* Ik D* m + a Ik D* I2l ) = ~^V d L ™, (38) 
and we find complete agreement for this term. 

The third term in (AFl) JI of Eq. (A. 8) of Ref.jHl] has the following correspondence in 



our notation 



One can prove that 



( AFl)(! d term -> Vjjhj sin /?ef (/) (39) 



hjDnMVn = V£t 
-j= smPi-VZXuUZ + X 2k U* 2 + tm6 w X lk U£) = V^ k c % (40) 

and we see that we again agree with BCRS. We notice here that our expression does not 
have the form factor (^(x, y, z). This form factor comes from the k 2 term in the integral 
where a loop with two fermions and one scalar is integrated. This part diverges and it is 
used to renormalize the parameters of the theory and could be safely ignored as we shall 
see when we compare with one of their e's later. 
The fourth term in (AFl) JI in their Eq. (A. 8) corresponds to our 

(AFi)S, term - Vjjhj sin (3ef ] (J) (41) 



One can prove that 



giV*^ - KjV^D^Vjj = -VJSS (42) 



and we find once again agreement. 

Next we compare with (AF|) J/ in Equation (A.9) of BCRS. The first term in (AF|) J/ 
has the following correspondence to our work 

term - -Vjjh* sin Pie'P (J) +e? 2) (J)) (43) 

The only difference that appear from the comparison is that the first term in BCRS should 
have an extra factor of e^ 3 . The second term in (AF^) JI of Eq.(A.9) in BCRS has the 
following correspondence to our work 

(AFl)£ d tcrm - -Vjjh* sin/3(e? 3) (J) + e? 4) (J)) (44) 
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One can prove that 



{®Ik D Ilj ~ llk D *I2j) = -^VJdN 



(PnDju + a* Jk D J2t ) = -j=V%jjr, (45) 
and we find agreement between our result and that of BCRS. 

The third term in (AF^) JI of their Eq. (A. 9) has the following correspondence to our 
work 

(Ai3)& tcrm - -VnhTj sin(3ef\j) (46) 



One can prove that 



g(U tl D nj - KjU i2 D I2j )Vjj = -V^ c 



DC 



-j= sin P(V2XlV a + X* 2k V l2 + tan 6 w X* k V i2 ) = -V^h, (47) 

and we find no difference between our equations and theirs for that term. The fourth 
term in (AF^) JI of their Eq. (A. 9) has the following correspondence to our work 



( AF 2 R )H term - -V n h\ sin Pe'P (J) (4* 



One can prove that 



gKjUaD^Vj! = V%§ (49) 

and comparing their result with ours we find no difference here either. To summarize, 
for the case with no CP phases we find complete agreement with the work of BCRS. 
However, for the case of CP- violation we have explicit gluino phase dependence. We note 
that the BCRS analysis did not take into account the CP violating effects in the Higgs 
sector. Specifically, it is now known that in the presence of complex phases in the soft 
SUSY-breaking sector, the three neutral Higgs mass eigenstates are mixtures of CP-even 
and CP-odd fields with production and decay properties different from those in the CP 
conserving scenarios. The vertices of these mass eigenstates are affected by this mixing 
moreover this mixing can lead to important effects in susy phenomena [2ZJ- We explain 
now in further detail exactly where misses these effects. Thus in , the assignment 
of the neutral Higgs as CP even h° and H° and CP odd A in Appendix A. 3 where 
different AF elements of Sec. 2 are calculated does not hold for the case with CP phases. 
Also the decomposition of the electroweak eigenstates in Eq. (3.8) of BCRS does not hold 



12 



for the CP violating case. Eqs. of sections (3.3), (6.1.2), (6.2) and Eq. (6.61) should 
also be modified to take this mixing into account. We note, however, that the CP mixing 
effects in the neutral Higgs sector do not affect the b — > 57 analysis in this paper. 

Continuing with the comparison of our work with that of BCRS we find that in Sec. 
(3) of j33], the authors did not take into account violations of isospin in their analysis 
as they are working in the approximation of Eq. (|13|) above. It is known that the effects 
of violations of isospin can be large, and if such violations were included, then in their 
Eqs. (3.3) and (3.14) the corrections A^Yd, A U Y^, A U Y U and A^Y^ should have a suffix i 
where % labels the element of an isopsin multiplet. In other words instead of the above 
four corrections one should have eight. Specifically the corrections that appear in Eqs. 
(3.35) and (3.41) are in general different from those in Eqs. (3.9) and (3.16). We also 
note that Eq. (3.37) of Ref. [3*4] is derived based on the assumption of isospin invariant 
loop corrections. 

However, compared to the analysis of the works of |H [H"J E2] the analysis of Ref. is 
more general as it takes into account more loops. Thus the analysis of Ref. [S3] specifically 
considered the case where two heavy fermions and one heavy scalar are running in the 
loops. So in their Fig. (9) they considered the additional important corrections to the 
charged Higgs boson couplings to quarks. The loop in Fig. (9a) of [B3j corresponds to 
our loop 2(vi) and the loop in Fig. (9b) corresponds to our loop 2(v). By looking at their 
expression for 5 a e' t (I) we note that there is no summation over the squark states. Thus 
the first term of 5 a e' t (I) corresponds to the case of squark j — 1 of our — ef (I). Using 
the following property of the form factor function H(x,y): 

m x m 2 m\ m\ m 1 mj mj 

— il \. — i—) = — H \— 1— ) l 5U J 

mjj 7713 7713 m,2 rri2 

one can make the comparison between our expression and theirs. So in the limit of 
vanishing Left-Right squark mixings, our expressions limits to that of [S3] . However there 
is a minor correction to their equation even in that limit. Thus in the first line of their 
expression for 5 a e' t (I), Zj should read Zj and their parameter a /j should be modified a 
little to be: 

a lj = Z*[s w Z% + c w Z 2 J} - V2Z u Z 3 Jc w (51) 
where the cos 9w factor in their last term is missing. The second term of 5 a e' t (I) correspond 
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to our — e t (/) with j — 2 squark. By noting that 

2 

7 tfc = -gX lk tand w (52) 

our expression reproduces their second term with the minor change that Zj should read 
Z 1 ^ and with the above new form of the parameter a /j of Eq. (|51jl . One can repeat the 
same analysis to compare 8},e' t (I) of [HI] with our — e' t (I)- Here also the analysis of 
Ref.[51j ignores squark mixing and their expression for 8^(1) corresponds to the squark 
of j = 1 in our equation. We need the fact that 

Pik = —{\s w Xi k - \c w X 2k ) (53) 

C\y 2 

to reproduce their expression with the minor change of the definition in Eq.([51|) above. 
We note here that the authors did not consider the part of the loop which has the form 
factor C2{x,y,z) as mentioned earlier. Finally, Eq. (5.6) of Ref.[3*4] calculates e' t (I) that 
corresponds to our — (e^ (I) +ef 3 \l)). Our expressions are more general and they limit 
to Eq.(5.6) of Ref.jHl] if we ignore in our formula, the first, third, fourth and fifth terms 
of e^\l) and by ignoring 19 terms in our e'} 3 \l) as well. We should notice here that they 
are using a different defintion of Ab. In this part of the paper they use Ab to be our — Ab as 
this could be seen from their footnote 3 of section 5.2. Also we note that their expression 
for e' t (I) is only valid for the CP conserving scenario. Thus for nonzero CP phases, fi in 
the first term should read fi* and Ab in the second term should read Al. However, the 
assignment of the Z^ matrix elements here is exactly like ours and is different from that 
in the DGG paper (See our note after our Eq.(|B3|)). 

Finally we should mention here that, apart from the small differences mentioned above 
with the approximate formulae of BCRS, our formulae should be rather considered as 
extensions not corrections of them. Next we compare our analysis to other earlier works 
where an effective Lagrangian similar to ours has been used. 

4.1 e' b (t) 

First we compare our analysis with the work of Demir and Olive (DO) [22] • We note that 
the e t b of DO is identical to our e' b {t). DO computed two one- loop contributions to e' b (t), 
which correspond to the contributions e^(i) and e^ 3 \t), in the limit of small squark 
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mixings and large tan (3. Our e' b (t) in this limit becomes 

2 2 r> m i m~ 

e b (1) (t) = -EE^I^flAHr^^,^) (54) 
f^ipi 3vr \m- g \ \m g \ 2 \m g \ 2 

which is the same as the first part of e tb in Ref . |32j . (We note that C q of DO is our D q 
and there is a typo in Ref.|32] m that their \Cf\ 2 should be \C? l \ 2 ). Next using 

k cot f3al k al k = h 2 t X; k Xi k (55) 



'(3) 



we find that our e b (t) in the limit assumed by DO takes the form 

1,2 a m "2 m~ 

To compare with Ref. [32] we define a t = h 2 /4ir and set C = X. One finds then that the 
overall sign of this term in DO is opposite to ours. As will be discussed later the overall 
sign of this term as computed in micrOMEGAs ^U] is also in disagreement with the sign 
given by DO, but in agreement with our sign as given above. Furthermore, in Ref. |32j 
(Co)4i(Co)3i should be (C^^Cq)^. Aside from these corrections, the results of Ref. [32] 
for e' b (t) for the parts computed, are in agreement with our result. 

Next we compare our results with the work of Degrassi et. al. (DGG) [3]. Eq.(15) 
of DGG can be obtained from Eqs. (|54I56|) of our analysis. To compare with the results 
of DGG we have to keep in mind that in the analysis of DGG, A b and m g are real. The 
relation between our X and the N of DGG is X T = N* and thus we note that their 
N4 a N* 3 should read N a ^N 3a . Moreover, one finds that the overall sign of the Yukawa 
contribution in DGG should be reversed to agree with our sign. 

The analysis of the micrOMEGAs group JU] takes into account all the six e' b (t) contri- 
butions but restricted to the case of real parameters and using certain approximations. 
The value of ejf^(t) and e^ 3 ^) is identical to the ones of DGG. The simplified formu- 
lae implemented in micrOMEGAs for e' b (t) was derived in Ref.j3H] and e^ 5 ^ 6 '* (t) was 
derived in Ref. |39j . 

We begin by displaying our e b 2 ^ (t) and e b ^ (t) in the limit of small squark mixings 



e 



2 2 n a* 1 m? 771- 

T\t) = EE^^ 3 lA 2 /|Aur^^(^,^) (57) 
~[~[^ix tan/3 \mg\ \m g \ 2 \m g \ 2 
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4 2 2 ii* h 2 1 777? 

4 W M - -EEE^Wli> B P«i^(^^) (58) 

fe=l 1=1 j=l ^ X fc Xfc Xfc 

To compare with the analysis of micrOMEGAs we keep in mind that in their work m g ,fi 
and A are all real. With this restriction our analysis is in agreement with Eq.(B.67) of 
Ref. [TOj specifically with the sign. (However, N* 4 in Eq.(B.67) should read N a4: ). Thus we 
support their disagreement with the work of DO and of DGG as stated after Eq.(B.67) 
in Ref 

Appropriately extended to the complex case the simplified formulae for (t) read 

, (5) = « Mfi (\Dnf H( W + W^W Wft m 

mi mi mi mi mi 



8sw7r V m? 




5TI//1 \ 1 1 Lr IILr llbT llir ULr 1 1 Or 

W ti ti (i *2 12 12 

and e'b 6 \t) read 

e;(6)(t) = a ^f ^M 2 | ^#(^, S-) + ^-H{^-, ^) | (60) 

These formulae are derived by using the corresponding formulae for and the decoupling 
limit. Moreover, one approximates the chargino masses by \x and M 2 and neglects mixing 
matrixes and U(l) contributions. We have checked numerically that they approximate 
the full formulae given in Eqs. (J81|84j) of Appendix B rather well over most of the complex 
parameter space. 



4.2 e bb 

In this section we carry out a similar analysis with the three works [H21 HI [HTj for the 
case of Em,. Comparing with the computation of DO we find that the QCD part given in 
Eq.(7) of Ref. [32] is the same as ours in the limit they are considering. To compare with 
the Yukawa part contribution we note that their Cl, Cr are related to our U and V as 
follows: C[ = V, and Cjj = U*. Then using 

g 2 — cot (3K t K b = ti (61) 

we find agreement with their analysis provide their (Cj^j is substituted by (Cjj)^- Next, 
comparing with the work of DGG, we agree with the QCD part of their Eq.(10) after 
taking account of the fact that they have no CP phases. To compare the contribution 
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of the chargino in their work with ours we note that their U is our U*. Also, V a2 in 
their work should read V* 2 . We note that there is also a disagreement between DGG 
and DO on this point taking into account that C L in Ref. [32] is V in Ref.[3] and C R in 
Ref.[32] corresponds to the matrix U of Ref.|3]. Finally we compare with the analysis 
of micrOMEGAs as given in Eq.(B.66) in Ref.|lUj. We agree with their result except that 
their V a2 should read V* 2 . The simplified formulae for 6$ + and + e$ extended 
to the complex case reads 

m , (2) _ 2a s (M SUSY ) (A b / tan f3 - fi*) < m\ 

e bb i* e bb — n M 12' 1 12/ V / 

67T Wig \ m g\ 

(3) ^ J4) _ yKM SVSY ) ^ /i/tan/3-4* rr\ mf 

66 



and 



e 



+ 4 } = E K 2 V ^'— ~^ RS ^t) (63) 

Xa X a Xa 

In micrOMEGAs the implementation of the terms 7 and 8 are given by e 7 bb = 2(e b 5 \t))* 
and e 8 bb = (e' fe (6) (t))* /2 using the results in Eas. (l59l(i()|) . 



4.3 e' t (s) 

First we compare our analysis with the result of DO, where et s corresponds to our e' t (s). 
DO only considered e t (s) and computed this in the limits mentioned in the preceding 
discussion. Our result in the same limits is given by 

3tt \m~ g \ \m g \ 2 \m g \ 2 

Using D s n ~ 1, D s \2 — 0, and m| = Q\ 2) we get exactly the e ts of Eq.(7) in Ref.[S2]- 
DGG only computed the tan/3 enhanced QCD and Yukawa terms; e' t (s) and e' t (s). 
Our e't (s) + e't \ s ) corresponds to their Eq.(16) reads 

2 2 9 i 2 /IT,? 

6; (1) ( S )+e; (3) ( S )=EE^-^|D iU | 2 |D^| 2 J- if( m * * ^ 



i^i 3vr |m § | |mg| 2 |m 



si 

/i 2 4 A* m 2 -. m f. 

~ : ^ s o 2 / — — A 3fc X 4fc |.D S 2i| -Dtij — o^' — 27 (^) 
fc=i *k xi xi 

Using the definition of their Eq.(17) there is a disagreement with the sign of the second 
part of their Eq.(16). Also, their Nl a N a3 should read N* 4 N* 3 . The micrOMEGAs group 
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only computes e t (s) and e' t (s) in our notation. Our approximation of the sum of these 
quantities gives 

+ 4 (2) (*) = E E ^V«V + r,)l^l 2 IA 23 fn ff (A> A) (66) 

~[~[ Sir tanp |m§| |^g| l m sl 

The analysis of the micrOMEGAs group does not include CP phases. Setting the phases 
to zero; /i* = /i, and £3 = etc. in Eq. (|66J) and taking into account that they use the 
opposite sign convention for e' t (s), we find that our result is in complete agreement with 
Eq.(B.68) in Ref.pJ]]. 

5 Numerical analysis and size estimates 

We now present a numerical analysis of our analytic results and also give a comparison 
with the previous works. In the following we compare three different methods for the 
calculation of the branching ratio of b — > 57, by using different computations of the e's; 

1. F: This is the full calculation of this work. 

2. S\. Here the e's are calculated using the simplified formulae found in the micrOMEGAs 
manual. Thus, we use the simplified formulae derived in Refs. jH ESI EHj appropri- 
ately extended to the complex derived in Sec. 4. Moreover, we correct the 
neutralino mixings terms entering in (t) and e' b (t) as stated in Sec. 4. 

3. S2: Here the e's are calculated using the simplified formulae of Ref . |32j . But with 
the corrections stated in this paper. 

In the previous section we exhibited the equivalence of our analysis and that of [SI] for 
the case with no CP phases, and showed that for the case with CP phases our analysis is 
more complete. Thus we use our analysis in the numerical computation since it is valid 
with and without CP phases. In the numerical analysis we take the SUSY scale to be 
the average of the stop-masses. We calculate the difference in percent to the full b — > S7 
calculation via the relation 

BR(b -> S1 ) F - BR(b -> S1 ) s 

BR(6^s 7 ) F 1 ' 

where S — S\^S%. In our numerical analysis we investigate several different supersymme- 
try breaking scenarios. These are 
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1. mSUGRA with real soft terms and complex SUGRA with universal value for the 
absolute soft terms. 

2. MSSM with real and complex soft breaking sector. 

In the analysis we scan over the parameter-space in order to find the qualitatively dif- 
ference among the schemes above and search for the parameter space where the allowed 
points satisfy the experimentally measured rate for b — > sj within 2 a, thus requiring 

2.3 x KT 4 < BR(6 -> s 7 ) F < 4.7 x KT 4 (68) 

In addition we check that all bounds on sparticle masses are satisfied, where we use the 
bounds given in Ref . |4Uj . Furthermore, we require the Higgs mass to be heavier than 110 
GeV in the real the theoretical error in the calculation of its mass is of order a 

few GeV. In the complex case we choose the lower bound 100 GeV for the lightest Higgs 
mass, as in this case there is a possibility for such a low mass being consistent with the 
LEP data jlj. This choice has little influence on our results. For the computation of the 
Higgs mass we use CPsuperH |4^] . 

Clearly, some of the contributions in the e's are numerically insignificant. We find 
that in e&& the contributions and ej;^ are small. Also the contributions e^ 3 ^ 4 ^) can 
be safely neglected, as the terms that would have been dominating are suppressed by 
the strange-quark Yukawa coupling. However, the contribution e'^'^^s), which has not 
been included in Si and S2 calculations of the rate for b — > s 7 gives sizeable contribution, 
capable of changing the rate by a few percent. The CKM elements V qq i that enter the 
analysis above are radiatively corrected and are calculated following the work of [3T]. 
Numerically the radiative corrections are found to be small in the part of parameter 
space investigated but the corrections could be significant in other parts of the parameter 
space. 

Before proceeding further we exhibit the dependence of the e's and the b — > 57 branch- 
ing ratio on phases. This is done in Fig. (5) which shows sharp dependence of these 
quantities on the phases. Specifically the analysis ob b — > 57 in Fig. (5) shows that the 
effect of phases can move the branching ratio for a given point in the parameter space 
from the experemintally forbidden area into the allowed corridor of values in Eq.([68j). 

In the following we discuss the different scenarios in detail. 
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5.1 mSUGRA and complex SUGRA with universalities 

In this section we carry out an analysis in the framework of extended mSUGRA model 
whose soft breaking sector is described by the parameters m Q ,A = \A Q \e ia , tan/3, /xo = 
|/io|e l6lM , rhi = \mi |e 1 ^ (i=l,2,3) where m is the universal scalar mass, mi is the universal 
gaugino mass, Ao is the universal trilinear coupling, and /io is the Higgs mixing parameter, 
while a, are the phases, all taken at the GUT scale. 

In the real mSUGRA case the scan is done by randomly selecting points within 
the following parameter-space; m G [200, 1000] GeV, my 2 G [200, 1000] GeV, A G 
[— 3mo, 3mo], tan/3 G [5, 55] and both signs of the \x parameter. For the complex case we 
also vary the five phases 9^, a a, £i, £,2, £3 within the range zero to 71. 

Our results are shown in FigJSJ We find a significant correlation of the increase of 
the differences with tan/3. This is very natural as the physical important parameter 
is etan/3 compared to one. Thus, in order for the e parameters to have a substantial 
influence tan/3 must be large. We see that in both cases the micrOMEGAs approximation 
is better (once we include the appropriate phases on their expressions) for the real case 
the differences remain below 2% using method S\ and for method S2 it is less than about 
4%. In the complex case while the Si approximation remain below 4%, the differences 
in the £2 approximation can reach 8%. This can be attributed to the fact that S2 does 
not include the electroweak contributions to e bb and e b (t). These contributions can induce 
a relatively large error at small values of mi/2 (indeed all the points with S2 ~ 8% 
correspond to m < 400 GeV and m 1( / 2 < 250 GeV). These results can be applied also to 
the super symmetric corrections to the b-quark mass, Sm b , which is given by e&&tan/3. We 
find that in the mSUGRA and complex SUGRA cases the simplification of S2 provides 
an accuracy of about 40% and the simplification of Si an accuracy of about 5%. 

We would like to stress the importance of using the correct signs and complex- 
conjugates. In Fig0we compare again the methods Si and 5*2 against our full calculations, 
but this time we use the original formulae, as presented in Ref . |32j and Ref.^Oj 4 - The 
difference is seen to increase substantial being as much as 15%. This can be understood 
since often there are cancellations among the epsilons contribution to the Wilson Coef- 
ficients. For instance the SUSY corrections to the W-contribution scale with e* bb — e' b (t) 

4 We note that the calculation in micrOMEGAs , is indeed performed correctly, despite of the error in 
the analytical formula. This is due to the fact that they use the a real N and thus allow for negative 
mass-eigenvalues, which is numerically a valid procedure in case of real parameters. 
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and this factor is much smaller than the e's themselves due to cancellation. Thus, having 
a wrong sign on one of the terms can cause a large effect on the SUSY correction. 

To compare the accuracy of the various methods of evaluation of the e's we focus on 
the point: 

m = 350, mi/2 = 200, A = 700, & = 0, £ 2 = 0.75, £ 3 = 0.5, 6^ = 0.6 a = tt , 

(69) 

where all the phases are given in radians and masses in GeV. Our results are shown 
in Fig. |H| For method Si the differences can be attributed to the simplification of the 
calculations of some terms and to neglecting the terms e t (s) and e t (s). We note 
that the inclusion of the electroweak contributions in e^ and e' b (t), is an improvement as 
compared to the simplified S2 method. 

Overall the corrections given in this paper to the BR(6 — * sj) are relatively small in 
the SUGRA scenario. Clearly, the SUGRA scenario constrains the MSSM mass spectra 
to have certain hierarchies. Moreover, due to the RGE evolution in SUGRA one normally 
finds the low-energy trilinear top term to be A t ~ — M 3 , unless the GUT scale A is 
very large. As we now discuss, these constraints in the SUGRA scenario give rise to 
various cancellations. The phase of the LO chargino contribution (see Eq. (J32)l ). assuming 
the hierarchy imposed in SUGRA scenarios, is given by Arg(/^4 t ) and the phase of the 
NLO chargino contribution is Arg(— //M3). The LO Higgsino contribution as well as W 
contribution are always positive. Thus, for the Higgsino and the chargino contribution 
to cancel against each other one needs Arg(fiA t ) ~ it. If such a cancellation occurs the 
SUSY corrections are allowed to be large. As noted in Sec. 3 all SUSY corrections scales 
with 1/(1 + e£ 6 tan/3) and thus these corrections will be large if e£ 6 tan/3 is close to minus 
one. The leading SUSY QCD contribution to e&& has a phase of Arg(^Ma). And this term 
is positive in the case of a negative chargino contribution. Thus, in mSUGRA one can 
never have a cancellation between the Chargino and the Higgsino contributions and at the 
same time have a negative value of e^. This is the main reason that the differences in the 
SUGRA case are rather small. Another reason is that in general in the SUGRA scenario 
the different contributions to the e's cancels against each other. Thus, for instance for 
the €bb correction, the leading SUSY QCD correction has the opposite sign as compared 
to the Yukawa and the electroweak contributions. Again this cancellation arises due to 
the relation between the trilinear top term and the gluino mass. Thus, in the SUGRA 
scenario the e's are numerically smaller than in the general MSSM scenario. 
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5.2 MSSM with real and complex soft breaking sector 

In the MSSM case the scan is done by randomly choosing the soft-masses in the range 
200 to 1000 GeV and the trilinear terms between —3000, 3000. We also run with plus and 
minus sign on the //-term, the trilinear terms and the gaugino masses. In the complex 
case we take the phases between to ir. We take the first and the second generation 
squarks and sleptons to be degenerate and we parameterize their masses by m s i. In the 
CP-violating case we use fairly heavy first and second generation, m s i = 2000 GeV, in 
order not to generate large EDM's. The first and second generations do not influence 
our calculation a lot, but they do enter in the evaluation of the chargino contribution to 
the Wilson coefficients. They also enter in the evaluation of e' t (s), but only in the terms 
e'} 3 \s) and (s), that can be safely neglected as they are numerically very small. 

In the MSSM case FigE] shows that the difference using method Si can be as large as 
60% and for method 5*2 we find roughly the same upper bound but the average difference 
is larger. Although on average the simplified formulae do a good job, there are cases 
with large errors in the rate for b — > 37. The large difference occurs in the MSSM case, 
particularly when there is a large sbottom mixing and/or large stop mixing. It is not 
difficult to realize this by looking at the formulae for the e's. Looking at the simplified 
formula for , one notice that these neglect the sbottom mixing, as compared to the full 
formulae. However, in the limit that the sbottom masses are equal the dominant term 
in the full formulae is invariant under sbottom mixing. Thus, in order for this correction 
to be important one needs a large sbottom mixing and a large sbottom mass difference, 
which occurs 'rarely'. As an example, in the formula for (t), the first term (neglected 
in previous works) is 

2 2 9 1 Til? TYl~ 

37r m b \m~ g \ \mg\ l \m g \ z 

The factor mt/mj, easily overcomes the suppression by cot (3. However, this term is zero in 
the limit of no sbottom or stop mixing. But, when the sbottom and stop mixing is large 
it can contribute significantly. Even the term with mf/m& cot fiD^D^ can give non- 
negligible contribution. For point (i) defined in Tablefl]we have the particular situation 
of both large stop as well as large sbottom mixing. We show the values of the e's and the 
rate of b — > 57 in Fig. lfTUj) . It is seen that in particular the value of e' b {t) is deviating from 
the value of the simplified formulae. Notice that this point with \x negative is excluded 
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with the simplified formulae, but allowed with the full calculation. The opposite might 
also occur as shown in Fig.([llj). 

The phase of the chargino contribution in the MSSM case depends on the mass hier- 
archies of the sparticles, and is thus less restricted than in the SUGRA scenario. Also, 
in the MSSM we are no longer confined to have a ~ n — £ 3 . Therefore, it is possible to 
have points where the chargino and the Higgsino contribution have opposite signs and at 
the same time have e^tan/? negative. Indeed all points with a difference of more than 
20%, assuming real parameters, have a value of e^tan (3 less than minus one half. For the 
point (i) plotted in FiglTUl the Chargino and Higgsino contribution are even larger than 
the SM contribution but as they cancels against each other, one finds a BR(6 — > 57) in 
agreement with experiment. 



Point 


M+ 


n 


M x 


M 2 


M 3 


Mq 


Mfj 


M b 


M L 


M E 


A t 


A b 


A T 


(i) 


450 


-950 


200 


-300 


400 


950 


900 


700 


300 


300 


2500 


1000 





(ii) 


300 


700 


200 


300 


-700 


550 


600 


500 


300 


300 


-1000 









Table 1: Values of the parameters for point (i) and point (ii) in GeV. The value of tan/3 
is not fixed and in both cases we have used m s i = 500 GeV. 

6 Conclusion 

In this paper we have given a more complete analysis of the next-to-leading-order contribu- 
tions which are enhanced by tan/3 factors. Such corrections affect the Wilson coefficients 
Cj and Cg arising from the W, Higgs H ± , and chargino x exchange contributions. There 
are twenty supersymmetric one-loop diagrams that contribute to these corrections. Some 
of these loops have been computed in previous works. In this paper we have given an 
analytic analysis of the full set of these corrections which involves computations of the 
six diagrams of Fig^ FigEI and FigJ3] each and the two diagrams of Fig0] The analysis 
presented here also includes the full CP phase dependence allowed within the general soft 
breaking in MSSM. The new analytic results of this paper are contained in Appendices 
B,C and D. In Sec. 4 we gave a comparison of the current work with previous analyses. In 
Sec. 5 a numerical analysis of the corrections was given and the effect of corrections found 
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to be significant specifically when there are large sbottom and stop mixings in the general 
MSSM case. The vertex corrections derived in this paper are relevant for a variety of 
phenomena where sparticles enter in the loops or are directly produced in the laboratory, 
such as Higgs decay widths and lifetimes and for the supersymmetric corrections to the 
b-quark mass. Since the analysis presented here takes full account of the effect of CP 
phases on b — > s + 7, it should serve as an important tool for testing supersymmetric 
models. 
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Appendix A 

Here we list the numerical values of the coefficients h iy h iy a^, 6, that appear in Eq.(8) jHj. 

,626126 56281 3 1 

hi = , , — , ,-0.6494,-0.0380,-0.0186,-0.0057 

v 272277' 51730 7 14' ' ' ' ; 

,14 16 6 12 
a% = ( 2^'2^'2^'~23' ' 086,-0 30 '"°- 899 ' ° 56) 
hi = (-0.9135,0.0873,-0.0571,0.0209) 

k = (0.4086,-0.4230,-0.8994,0.1456) (71) 
Appendix B - Analysis of e' b (t) 

The analysis of e' b (t) as well as of e' t (s) and of e&& depends on the soft breaking pa- 
rameters. We shall carry out the analysis in the framework of MSSM which has a pair of 
Higgs doublets with Higgs mixing parameter fi which is in general complex, assuming a 
general set of soft breaking parameters. Specifically we will assume for the e analysis a 
general set of squark masses, and of trilinear couplings A q which we assume in general to 
be complex. Similarly, we assume the gaugino masses fhi (i=l,2,3) to be complex. Thus 
our analytic analysis will not be tied to any specific model of soft breaking. There are six 



24 



different loop diagrams that contribute to e' b {t) so that 



4(*) = E4 (i) (*) ( 7 2) 



1=1 

We exhibit now each of the above contributions. 
From Fig^i) we find 

2 «s if, r»* r m * 



4 J (*) = "EE -^e^D* b2j D tu [-* cot f3A t D blj D; 2i + ^D b2j D* tll + m t cot /3D b2j D t * 2i 



i=l j=l 



m 2 m 2 1 ^1 

+ _L cot ^ D * -^sMcos/^A'iJi — r#(Ht^T^k) ( 73 ) 

m fe m b \m g \ \m g \ z {rrigf 

where D q is the matrix that diagonalizes the squark mass-squared matrix, i.e., 

D\Mp q = &mg{Ml,Ml) (74) 

From Fig^ii) we find 

4 (2) (*) = r^EE ^e^DJyDtti^Ky^ + — // cot (3D blj D* m + m t D b2j D* t2i 

tan p i=1 . =1 07T 

m 2 1 m- 

+m b D blj D* m - — — cos 2 p^ Wj £>• J if (75) 

This diagram is not enhanced by tan/3. From Fig^iii) we find 

£?\t) = 2 j^jzjzi— cot PA t D bl] Dl 2l + ^D b2j Dl u 

k=li=lj=l m b 

2 2 

+m t cot (3D b2] D* t2l + ^ cot (3D bl] D* m - ^ sin /5 cos (3D blj D* tll ] 
m b m b 

« k D* blj - ll k D* b2j ){(3 tk D m + 3 L ) (76) 

x fc x° x£ 

In the above a, p\ and 7 for the b and t quarks are defined so that 

gm&(t)X 3 (4) fc 



b ^ k 2mw cos /3 (sin 0) 
Pb(t)k = eQb(t)Xi k H ^— X 2 fe(T 36 ( t ) - Qfe(t) sin 2 

COS C7VK 



- y' ^Qfe(t) sin2 v > (77 ^ 
7b(t)k - eQ b (t)A lk X 2k (77) 

COS U\y 
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where Q b{t) = -§(§) and T 36(t ) = -|(|) and where 

X' lk = X lk cos #w + X 2k sin # w , X' 2k = —X lk sin # w + X 2k cos W (78) 
and X diagonalizes the neutralino mass matrix. 

X T M x o X = diag (m x o , m x o , m x o , m x o ) (79) 
From Fig^iv), which is non-tan (3 enhanced, we find 



taI1 P fc=li=l i=l m & 



2 

-m t D b2j D* t2i + m b D blj D* tli - '-^ cos 2 (3D blj D 



m b 



tin 



H k D* blj - r bk D* b2j )(P tk D tli + a * tk D t2l )-^—H(^ ^) (80) 

107T 77l v o /7i_o ^,,0 

x fc x fc x k 



From Fig^v) we find 



4 2 2 

4 (5) (*) = E E E cos p^f^i,] (V2X^Va + X* 2k V l2 + tan 8 w X* k V l2 ) 

k=l i=l j=l mfo 

1 m-m x o m 2 m\ 

* (f3 tk D tlj + a * tk D t2j )H(-^, -f ) (81) 



16tt 2 mf tiJ v ? > 2 

In the above £7 and are the matrices that diagonalize the chargino mass matrix 

U*M X+ V~ l = diag(m Y +,m y +) (82 

and K b is given by 



K b = mb (83) 
y/Imw cos p 



Finally, from FigGJvi) we get 



4 2 2 



4 \t) = -E EE V2g—cosP(UaD blj - K b U i2 D b2j ) 
k =i i=i i=i m & 

(V2X 4 *^i + X* fc ^ 2 + tan6 w X* lk V l2 ) 
1 m-m x o m 2 m 2 
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Appendix C - Analysis of e' t (s) 



Next we look at the e' t (s) analysis. Here we have 

<(s)=j:4%) (85) 



The individual contributions e t (s) are exhibited below. 



f '(0 

i=i 

From Fig|2fi) we find 

2 2 r> 

<=i j=i 3?r m * 

777 2 777^ 1 777 ? 777 r 

+ ^tan/? J D sllJ D* - — sin {3 cos [3 D sli D* }—H(-^-, -^-) (86) 
Wlf J m t \ m g\ \ m g\ \ m g\ 

From FigO^ii), which is non-tan f3 enhanced, we find 



m 2 1 771? m ? 

+m s A^ + m t D slj ^. - sin 2 /JD. M I?y —Iff-^, r-%) (87) 

From FigEfiii) we find 

4 2 2 



e ; (3) f S ] 



777 

- 2 E E Et— ^/3A* s D s2l D* tlj + ^*D s1i D; 2j + m s t^^D m D* t2j 

k=n=ij=i m t 

2 2 

+^tan / 3L> sli L>* - ^ sin /3 cos ^ sli D* ] 
m 4 J m t 

1 1 m 2 m f 

(^:« + - itkD t2j )—— h(-j- , -f ) (88) 

lOTT /7£ y '''',,0 ,0 

Afc X k x fe 



From Fig|2fiv) we find the non-tan (3 enhanced contribution 



2 * * * „ _ m. 



4 (4) (^) = -^EEE^^ + ^tan^Ai, 
tan P fc=ii=i j=i m * 



2 

777 

+m s D s2i D* t2j + m t D sli D* - sin 2 [3D sli D* tlj \ 

1 1 m 2 m ? 
(oteDty - TtoAtfXfo^M + « sfeJ D: 2l )-— F(-f , ) (89) 

1D7T 777yO 777 777 
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From Fig|2fv) we find 



where 



From FigO^vi) we find 



4 2 2 



EEE^^sin/3^^] 

fc=i i=i j=i m * 

(-v^f/* + X 2k U* 2 + taa^X lfc C^) 



m 2 _ m 2 



1 m y om v - 

(/3* k D alj . + a Bfc D; 2j -) 1fi 2 Xfc *' #(^S -£) (90) 
7 I67H m 5 2 ml mi 



Kt = (91) 



V^miy sin /3 



4 (6) («) = " E E E V^tf— sin Ay - K t V* 2 D* t2j ) 

1 m.om- m 2 m 2 
c^Ay - 7«fcA 2j ) 1fi 2 % X ' (92) 



ij tj tj 



167T 2 Mr 771- 771 

Appendix D-Analysis of 



We proceed now to compute the e bb . It is given by 

^ = E4? (93) 

i=l 

We exhibit the individual contributions below. 
From FigEfi) we find 

e bb = -LL^ e ^ D bii D b2 3 ¥ D b\iD h2j + — 94 

37r J m b cos9 w 

11 1 1 ml m l 

{(-- + - sin 2 6 w )D m D* blj - - sin 2 6 w D b2l D* b2j } sin/3] 1 _2y 

2 3 J 3 J |m§| {mgl \ m g\ 

From FigOfii) we find the SUSY QCD non-tan j3 enhanced contribution 

1 ^ — > ^ — ^ 2cv ; 

37r 



-EE ^^« 2j [-4« 2j - m b {D* blj D m + D* b2 -D b2i } 
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m z m w cos/3 1 1 . 2 1 . 2 

7 — 1( — o + o sm Vw)D bl iD blj - - sin 6 w D m D b2 A cos/3 

m 6 costly 2 3 J 3 J 

p-r^p-Vp-^) (95) 
l m !?l l m <?l \ m g\ 



From FiglSfiii) we find 



4? = ~ E E E # 2 [- — cot P A t D t2iD*tij - — ^ P{D tll D* tlj + D t2i D* t2j } 
i=i j=i k=i m & mb 



m z m w cos/3 1 2 . 2 2 . 2 . 

~ m ;^ _ ^9~q sm ^w)AiiAij + o sm 9 w D t2l D t2j 

TTlb COSt7jy Z O J 6 J 



1 1 m~ m~ 

- K t V k \D^{K h U* k2 D tlj )—- 2 - 1 #(,— V ( 96 ) 

lo^ m~+ \m~+\ z \m~+r 



From FigOfiv), which is non-tan j3 enhanced, we find 

4? = -^|g|:/ 2 e-w- 



m z m w cos/3 1 2 . 2 2 . 2 

m b cosv w 2 3 7 3 7 

11 m? m~ 

(y fe *L^ - ^^(^t^,^) .^ (ur^, ur^) (97) 



|m^+ \m~+r m--, 



From FigOfv) we find 



i=l j = l fc=l 

1 1 

2 + 3 " 



m b cos 6*iy 



1 • 2 fl \ n n * 



{( _ o + o sin Qw)D h iiD* blj - - sin e w D b2i Dl 2j } sin/3] 



,2 /m? 



1 1 m| mg. 

(a bk D blj - -f bk D b2j )(p; k D* bli + a bk D* b2i ) ^— 2 , #( ,_ _' |2 , _ J I2 ) (9* 



m^o I ' im^o | 2 ' im^o | 



From FigOfvi), which is non-tan /3 enhanced, we find 



4? = AEEE2 - m 6 { + A«Z>Sy } 

tanp i=1 j=lk=l 

m z m w cos/3 1 1 . 2 1 . 2 

77— {(-t + 77 sm 9 w )D m D blj - - sm 9 w D m D b2 A cos/3 

m b costV 2 3 J 3 7 

1 1 ™i m i 

(o*Ay - T*D*j)Wrtu + a ^ D ^T^M H{ \^> M } (99) 
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From FigEJi) we find 



2 2 2 j _|_ 

4? = EEEV^ cot - Q*- cos (3 - E* ] (100) 
i=l j=i k=l " Lb *' n w 

I \m~+\\m~+\ |m-+| 2 \m-+\ 2 
KDU - K t V; A){ K b U h D tlt) — 2 ^-^H { ^, 

From FigE^ii) we find 



i=i j=i fe=i 



b2fc J 



1 l m x°ll m x?L. J m v°| 2 l m x?l 



+ ^k)^ A mg -^H (ioi) 

In the above Q, R, Q" and R" are defined by 

Qij = J^U i2 Vji 

R* = ^[^UaVji + iSUaVp] (102) 



and by 



gQy = -[X^gX^-g'X^ + ii^j)} 
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Figure 1: Set of diagrams contributing to e£(i) 
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Figure 2: Set of diagrams contributing to e' t (s) 
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Figure 5: Variation of the e's and BR(6 — > 37) with tan/? for the parameters m = 
400 GeV, mi/2 = 200 GeV, Aq = 700 GeV and all the phases set to zero (solid lines) and 
with the phases (in radiants) set to £1 = 0, £2 = 0.8, £3 = 0.9, 6^ = 0.6, «o = re (dash 
lines) or £1 = 1, £2 = 0, £3 = 0.3, 9 ^ = 1.2, «o = tc/2 (dot-dash lines). The thick lines 
represent real values while the thin lines correspond to imaginary parts. 
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Figure 6: The percentage difference between the approximate formulae of Si and S2 and 
our full calculation in the SUGRA scenario. The left graph is the case with no phases, 
and the right graph is the case with phases. Each group contains about 1800 models 
where each point in the parameter space defines a model. 
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Figure 7: The percentage difference between the non-corrected approximate formulae as 
they appear in Ref.[32] and Ref.fTU] and our full calculation in the real mSUGRA case 
which has no phases (left) and in the complex SUGRA case which has phases (right). 
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Figure 8: Variation of the e's with tan/3 for the parameters of Eq. (|69|) and the corre- 
sponding prediction of BR(6 — > 57) using S±, S 2 and the complete calculation (F). 
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Figure 9: The percentage difference between the approximate formulae of Si and S2 and 
our full calculation in the MSSM scenario. The data sets contain about 1000 models. 
We only plot points that are experimentally acceptable. Left graph is the case with no 
phases, and right graph is the case with phases. 
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